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Abstract 

We  report  the  first  direct  experimental  observation  of  the  ablative  Richtmyer- 
Meshkov  instability.  It  manifests  itself  in  oscillations  of  areal  mass  that  occur  during  the 
shock  transit  time,  which  are  caused  by  the  rocket  effect  or  dynamic  overpressure 
characteristic  of  interaction  between  the  laser  absorption  zone  and  the  ablation  front. 
With  the  4  ns  long  Nike  KrF  laser  pulse  and  our  novel  diagnostic  technique 
(monochromatic  x-ray  imaging  coupled  to  a  streak  camera)  we  were  able  to  register  a 
peak  and  a  valley  of  the  areal  mass  variation  before  the  observed  onset  of  the  Rayleigh- 
Taylor  growth. 

PACS  numbers:  52.57.Fg,  52.70.La,  47.20.-k 
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Laser  target  distortion  due  to  the  Rayleigh-Taylor  (RT)  instability  remains  a 
critical  issue  of  the  inertial  confinement  fusion  (ICF)  program  [1,  2],  The  perturbations 
amplified  by  the  RT  instability  grow  from  the  so-called  seeds  which  are  linearly  related 
to  the  non-uniformities  of  the  target  and  the  laser  beam  or  the  x-ray  radiation  field  driving 
it.  The  growth  rates  of  the  RT  instability  have  been  measured  in  numerous  experiments 
starting  from  the  mid-1980s  (e.  g.,  see  Refs.  3).  Much  fewer  experimental  results  are 
available  to  quantify  the  linear  relationships  between  the  measurable  initial  non¬ 
uniformities  and  the  seed  mode  amplitudes.  A  good  understanding  of  these  relationships 
is  necessary  to  ensure  a  reliable  simulation  of  the  overall  perturbation  growth,  and  hence, 
of  the  target  performance. 

Most  of  our  knowledge  in  this  field  still  comes  from  simulations  and  theory.  The 
seed  amplitudes  are  known  to  form  during  the  early-time  period  that  includes  a  shock 
wave  transit  from  the  front  to  the  rear  surface  of  the  target,  and  a  rarefaction  wave  transit 
in  the  opposite  direction.  During  this  time  interval,  the  areal  mass  perturbations  caused  by 
all  sources  of  non-uniformity  are  expected  to  oscillate  [4-9].  The  physical  mechanisms 
driving  the  oscillations  depend  on  whether  the  perturbations  are  initially  at  the  front 
surface  of  the  target  (laser  imprint,  front  surface  roughness)  or  at  its  rear  surface  (feedout, 
see  [9]).  Here,  we  limit  ourselves  to  the  former  case,  where  the  oscillations  are  caused  by 
the  rocket  effect,  or  the  dynamic  overpressure  [6,  8,  10,  11].  The  oscillatory  behavior  is 
consistent  with  the  expression  for  the  growth  rate  r  of  ablative  RT  instability,  which  has 
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Here  g  is  the  acceleration,  rD  « 1  is  the  effective  blowoff-plasma-to-ablation-front 
density  ratio,  estimated  for  a  given  perturbation  wavelength  A  =  2n  /  k  by  the  formulas 
given  in  [8,  11],  A  =  (l  -  rD )/ (l  +  rD  )  is  the  corresponding  effective  Atwood  number, 
va  =ml  ps  is  the  ablation  velocity  ( rii  is  the  mass  ablation  rate,  ps  is  the  shock- 
compressed  target  density).  The  main  difference  between  (1)  and  all  the  versions  of  the 
Bodner-Takabe  formula  (see  [12]  and  the  references  in  [1-3])  is  the  presence  of  the 
negative  term  under  the  square  root.  The  so-called  rocket  effect  described  by  this  term 
emerges  because,  as  found  in  Ref.  12,  the  ablation  front  is  an  isotherm.  When  it  is 
perturbed,  and  a  part  of  it  gets  closer  to  the  hot  laser  absorption  zone,  the  temperature  at 
the  ablation  front  does  not  increase,  but  the  temperature  gradient  in  its  vicinity,  VT  , 
does.  This,  in  turn,  increases  the  local  heat  flux  to  the  ablation  front,  -kVT  ,  hence,  the 
rate  of  mass  ablation  from  it,  thereby  increasing  the  ablative  pressure  and  pushing  this 
part  of  the  ablation  front  back.  The  physics  of  this  rocket  effect  is  explained  in  detail  in 
Refs.  6,8,10.  The  rocket  effect  rather  than  the  mass  flow  through  the  ablation  front 
[described  by  the  last  term  in  the  right-hand  side  of  (1),  see  Ref.  12]  determines  the  cut¬ 
off  wavelength  of  ablative  RT  instability  at  high  Froude  numbers.  Formula  (1)  giving  the 
rate  of  exponential  growth  [10,1 1]  is  formally  applicable  even  when  the  acceleration  is 
low,  or  the  wavelength  is  short,  and  such  a  growth  does  not  take  place  [6,  7].  If  we 
substitute  the  value  of  g  =  0 ,  appropriate  for  the  shock  transit  time,  into  (1),  then  the 
contribution  of  the  square  root  term  is  purely  imaginary,  T  =  ?  where 

Q.  =  Akva  /  r^2  (2) 

is  the  oscillation  frequency  [6] .  These  oscillations  during  the  shock  transit  time  were  first 
predicted  in  our  simulations  [5],  then  explained  by  the  theory  of  Ref.  6  (where  the  term 
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ablative  RM  instability  was  introduced),  and  later  confirmed  by  many  other  simulations 
[7,  8],  but  never  actually  observed  experimentally. 

One  of  the  problems  associated  with  observation  of  the  ablative  RM  oscillations  is 
their  relatively  low  frequency.  This  is  why  these  oscillations  have  not  been  detected  by 
the  face-on  measurements  in  earlier  experiments  [13].  For  a  30  pm  perturbation 
wavelength  in  a  plastic  target  irradiated  by  high  laser  intensity,  ~80  TW/cm  ,  va  =  1  to  1.5 
pm/ns,  rD  =  0.04  to  0. 1  [8,  11],  Eq.  (2)  predicts  an  oscillation  half-period  between  2  and  4 
ns.  The  driving  laser  pulse  has  to  be  even  longer  than  that,  if  we  want  to  observe  the  RT 
growth  that  follows  the  early-time  RM  phase.  This  would  provide  the  data  to  validate  our 
codes  that  model  the  early-time  RT  seed  generation.  To  ensure  a  sufficiently  long  shock 
transit  time,  the  plastic  targets  must  be  at  least  40  to  60  pm  thick.  Observation  of  the 
oscillatory  behavior  of  mass  variation  in  such  targets  requires  direct,  high  spatial 
resolution  face-on  measurements,  registering  relatively  small  changes. 

We  have  performed  experiments  in  this  regime  to  observe  the  evolution  of 
ablative  RM  instability  in  laser-driven  targets;  we  present  here  a  subset  of  our  data  to 
summarize  our  results.  Our  experiments  were  performed  with  the  Nike  KrF  laser  [14] 

(  =  248  nm).  The  laser  pulse  (~1400  J  in  37  overlapping  beams)  was  focused  to  a  spot 

750  pm  FWHM  in  diameter,  with  a  flat  central  region  400  pm  in  diameter,  producing  an 
intensity  up  to  ~8xl0  W/cm  .  In  this  set  of  the  experiments  the  main  Nike  was  4  ns 
long,  satisfying  the  previous  pulse  duration  requirement.  The  time-averaged  rms  spatial 
variation  for  overlapped  37  beams  has  been  measured  to  be  less  than  0.2  %;  hence  the 
target  is  not  severely  distorted  by  the  laser  imprint  during  the  shock  transit  time. 
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The  diagnostic  setup  for  the  RM  measurements  (Fig.  T)  is  a  modification  of  the 
Nike  monochromatic  x-ray  imaging  system  based  on  Bragg  reflection  from  spherically 
curved  crystals  [15].  An  addition  of  a  streak  camera  to  this  system  made  it  possible  to 
analyze  continuous  time  behavior  of  mass  variation,  which  is  necessary  to  reveal  the  non¬ 
monotonic  evolution  of  the  processes  under  study.  The  energy  of  12  Nike  beams,  ~500  J, 
is  delivered  to  a  silicon  backlighter  target,  producing  x-rays  that  backlight  the  main  target 
for  about  5  ns.  The  spherically  curved  quartz  crystal  selects  the  resonance  line  of  the  He- 
like  Si  (1.86  keV)  and  projects  a  monochromatic  image  of  the  target  on  the  slit  of  an  x- 
ray  streak  camera.  We  used  quartz  crystal  with  the  cut  1011  and  radius  of  curvature  200 
mm.  With  magnification  of  20  it  makes  the  crystal-to-detector  distance  about  2  m.  Our 
system  was  able  to  position  the  target  image  on  the  entrance  slit  of  the  streak  camera  with 
an  accuracy  that  corresponds  to  30  pm  on  the  target.  This  setup  provides  spatial 
resolution  in  one  relevant  direction,  along  the  wave  vector  of  the  a0  sin  kx  ripple  on  the 
target  surface,  producing  a  streak  record  like  that  shown  in  Fig.  1(b). 

Spatial  resolution  of  the  x-ray  optical  system  was  tested  with  the  help  of  x-ray 
film  as  a  detector  and  was  found  to  be  of  order  6-7  pin.  High  throughput  of  the  x-ray 
optical  system  allows  us  to  magnify  the  images  by  a  factor  of  20  to  compensate  for  the 
modest  resolution  of  the  streak  camera.  Despite  the  high  magnification  of  the  system,  the 
overall  spatial  resolution  was  still  limited  by  the  streak  camera  and  was  estimated  by  test 
pictures  of  the  composite  mesh;  see  Fig.  1(a).  The  MTF  of  the  entire  diagnostic  system 
was  obtained  by  imaging  a  knife-edge  target  and  verified  by  imaging  an  undriven  target 
with  predetermined  amplitude.  The  MTF  at  A  =  30  pm  and  45  pm  is  0.4  and  0.7, 


5 


respectively.  This  test  also  gives  us  an  experimental  error  bar  for  the  mode  amplitude 
measurement,  which  does  not  exceed  ±  lpm  in  most  cases. 

We  used  40  pm  to  90  pm  thick  CH  targets  rippled  on  the  front  side  with 
perturbation  wavelength  A  =  30  pm.  This  was  short  enough  to  make  the  oscillation  half¬ 
period  fit  into  the  Nike  pulse  duration,  but  at  the  same  time  was  resolvable  by  our 
diagnostics.  The  initial  ripple  amplitudes  a0  were  approximately  either  0.5  or  1  pm.  This 

combination  of  perturbation  wavelength  and  initial  amplitude  might  not  be  adequate  for 
an  experiment  aimed  at  measuring  the  linear  RT  growth:  with  ka<>  =  0.1  to  0.2,  we  are  not 
far  from  the  non-linear  stage  from  the  very  beginning.  Here,  however,  we  are  primarily 
interested  in  the  early-time  ablative  RM  phase,  when  the  perturbations  of  all  wavelengths 
oscillate  rather  than  grow  exponentially,  so  that  the  higher  harmonics,  being  initially 
small,  remain  small  until  the  onset  of  the  exponential  RT  growth.  The  streak  records  were 
taken  with  time  resolution  of  170  ps,  sufficient  for  the  characteristic  times  of  interest, 
which  are  no  less  than  0.5  ns.  The  large  flat  top  (400  pm)  of  the  laser  focal  spot  observed 
with  the  large  field  of  view  (500  pn)  gives  us  more  ripples  available  for  Fourier 
transform  analysis,  thus  ensuring  confidence  in  determining  the  dominant  mode  s 
amplitude.  We  present  the  Fourier  amplitude  of  mass  variation  [pm_g/cm3]  in  microns, 
referring  it  to  the  density  of  solid  plastic,  1.07  g/cm3.  This  amplitude  should  not  be 
confused  with  displacement  of  the  ablation  front.  The  observed  mass  variation  also 
includes  contributions  from  the  rippled  shock  front  and  the  density  non-uniformity 
between  the  shock  and  ablation  fronts.  Our  face-on  diagnostics  does  not  allow  us  to 
evaluate  any  of  these  contributions  separately.  Theory  and  simulations  indicate  that  the 
contribution  of  the  rippled  ablation  front  rapidly  becomes  dominant  [5-8], 
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Fig.2  shows  an  original  streak  record  (top)  and  the  time  dependence  of  its  MTF- 
corrected  Fourier  component  (bottom,  peak-to-valley  amplitude)  taken  from  a  47  pm 
thick  rippled  target,  a0  =  0.5  pm,  irradiated  with  peak  intensity  of  5x1 013  W/cnr,  no  low 
energy  foot.  The  thickness  of  the  shaded  area  approximately  corresponds  to  the 
experimental  uncertainty.  The  peaks  and  valleys  of  mass  variation  are  clearly  visible  on 
the  streak  record.  Both  light  and  dark  stripes  are  less  discernible  between  2.5  and  3  ns,  as 
if  they  were  smeared  in  a  vertical  band  ~  0.5  ns  wide,  indicating  a  minimum  in  the 
Fourier  amplitude.  The  fast  rise  in  mass  variation  following  the  RM-related  minimum  is 
the  observed  RT  perturbation  growth.  Note  that  the  late-time  peak  at  about  4  ns  occurs 
when  the  laser  pulse  is  over  and  the  target  decompresses. 

Figure  3  compares  the  results  for  the  above  shot  with  a  thicker  target  (63  pm 
thick,  a0  =  1  pm)  at  about  the  same  peak  intensity.  The  shock  transit  takes  longer  for  the 
thicker  target,  so  the  RM  phase  is  longer  by  about  a  factor  of  1.5.  Hence,  the  RT  growth 
starts  later,  and  is  slower  due  to  the  higher  mass  (lower  acceleration)  of  the  target.  By  the 
end  of  the  laser  pulse,  the  mass  variation  amplitude  grows  to  approximately  the  same 
value  as  in  the  thinner  target,  despite  having  started  from  twice  the  initial  value. 

Figure  4  shows  the  data  obtained  for  the  thickest  target  that  we  probed  with  our 
diagnostics,  90  pm,  2«0  =  2  pm.  In  this  shot,  the  laser  beam  had  a  3  -ns  long,  1%  foot  at  t 
<  0,  which  caused  some  early  growth  of  mass  variation  (not  shown  in  the  figure)  due  to 
the  rippled  shock  wave  and  to  the  perturbation  of  the  ablation  front  to  a  ~4  pm  amplitude. 
Here,  the  shock  transit  takes  still  longer,  and  the  acceleration,  when  it  starts,  is  lower, 
hence  only  the  onset  of  the  RT  growth  is  seen  (the  laser  pulse  ends  at  4  ns).  Still,  the 
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growth  of  the  dominant  mode  to  the  peak  and  its  subsequent  decay  before  the  onset  of  RT 
are  clearly  observed.  This  result  is  compared  to  a  simulation  performed  in  two 
dimensions  (2D)  using  the  fast2d  hydrocode  developed  at  the  Naval  Research 
Laboratory  [16]  (more  details  about  the  code  and  further  references  are  given  in  Refs.  5, 
7).  Qualitatively,  the  simulation  predicts  exactly  what  we  observe  in  all  shots. 
Quantitative  agreement  is  also  quite  reasonable. 

In  summary,  we  have  observed  for  the  first  time  a  non-monotonic  evolution  of 
areal  mass  perturbations  in  an  ablatively  driven  target  during  the  shock-rarefaction 
transit.  The  perturbations  experience  a  half-oscillation  due  to  so-called  ablative  RM 
instability,  followed  by  the  onset  of  exponential  RT  growth.  The  initial  perturbations 
were  single-mode  ripples  on  the  front  surface  of  planar  targets,  corresponding  to  a  mode 
of  the  outer  surface  roughness  of  a  spherical  direct-drive  target.  The  physics  of  the  RM 
oscillations  excited  by  either  surface  ripples  or  single-mode,  constant  phase  laser  beam 
non-uniformity,  is  identical  [5-8].  Using  the  former  case  to  validate  the  codes,  the 
capacity  for  modeling  the  latter  can  be  improved. 

Of  course,  it  would  be  of  great  practical  importance  to  observe  the  ablative  RM 
oscillations  in  a  solid  deuterium  target.  Our  present  diagnostics  cannot  see  mass  variation 
in  deuterium,  which  does  not  absorb  keV  x-ray  quanta.  There  is,  however,  an  attractive 
opportunity  of  continuing  this  line  of  experiments  with  deuterium-wicked  plastic  foam 
targets,  as  in  Ref.  17.  The  oscillation  frequency  in  deuterium  should  be  higher  than  in 
solid  plastic  due  to  its  lower  density  [higher  va  in  (2)],  making  it  possible  to  observe  more 
than  one  half-cycle  [5,  6,  8],  whereas  the  total  mass  variation  would  be  detectable  via  the 
observed  changes  in  the  areal  density  of  the  foam  component. 
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Figure  captions 


Fig.  1.  General  scheme  of  the  experiment,  (a)  Test  images  of  a  150-to-400  lpi  composite 
mesh  obtained  on  x-ray  film  (quadrants)  and  on  an  open  photocathode  of  the  streak 
camera  (insert  in  the  middle),  (b)  An  actual  streak  record  showing  the  RT  growth  of  areal 
mass  variation  —  peaks  and  valleys  are  more  pronounced  at  late  time. 

Fig.  2.  Streak  record  (top)  and  the  peak-to-valley  amplitude  of  the  dominant  Fourier 
mode  A  =  30  pm  of  the  areal  mass  variation  (bottom)  for  a  47  pm  thick  target  (peak-  to- 
valley  amplitude  2 a0  =  lpm).  The  minimum  due  to  the  ablative  RM  oscillation  is  clearly 
visible  on  both  images. 

Fig.  3.  Time  histories  of  the  areal  mass  variation  (peak-to-valley)  compared  for  47  pm 
(2(7o  =  1  pm,  A  =  30  pm)  and  a  63  pm  (2«0  =  2.0  pm,  A  =  30  pm)  thick  targets  irradiated 
at  5xl013  W/cm2  (solid  and  dotted  lines,  respectively).  The  oscillation  is  followed  by  the 
RT  growth,  which  peaks  when  the  laser  turns  off  and  the  target  starts  to  decompress. 

Fig.  4.  Time  histories  of  areal  mass  variation  (peak-to-valley)  for  a  90pm  thick  target 
(2(7o  =  2  pm,  A  =  30  pm)  at  4.5xl0lj  W/cm2 —  the  RM  oscillation  is  still  observable,  but 
only  onset  of  the  RT  growth  is  seen.  This  specific  shot  was  taken  with  3  ns,  1%  foot.  The 
thickness  of  the  shaded  area  corresponds  to  the  experimental  uncertainty  of  the  mode 
amplitude  measurement.  Experiment  is  compared  to  a  simulation  performed  with  the  help 
of  FAST2D  hydrocode. 
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